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Abstract 

The possibility to observe the nuclear modification of the gluon distribution at small- 
X (gluon shadowing) using high-p^ prompt photon production at RHIC and at LHC 
is discussed. The per-nucleon ratio, a{p + A — > 7 + X) /(vl x a{p + p — > 7 + X )), is 
computed for both inclusive and isolated prompt photons in perturbative QCD at 
NLO using different parametrizations of nuclear parton densities, in order to assess 
the visibility of the shadowing signal. The production of isolated photons turns out 
to be a promising channel which allows for a reliable extraction of the gluon density, 

A A 

R^, and the structure function, R^ , in a nucleus over that in a proton. Moreover, 

G F2 

the production ratio of prompt photons at forward-over-backward rapidity in p~K 

A A 

collisions provides an estimate of R^ (at small x) over R^ (at large x), without the 
need of p-p reference data at the same energy. 
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1 Why and how to probe gluon shadowing 



Gluon distributions in a proton and in a nucleus are fundamental ingredi- 
ents in order to compute, within perturbative QCD (pQCD), hard-process 
observables in proton-proton, proton-nucleus, and nucleus-nucleus collisions. 
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Thanks to years of detailed experimental studies at HERA and Tevatron sup- 
plemented by important theoretical developments in global fit analyses, the 
gluon distribution G''{x,Q'^) is fairly well known, say within a few percent 
accuracy, in the range x ~ 10~^-10~^ and ~ 10-10^ GeV^ [1] which is 
precisely the kinematical domain covered by most hard processes at the LHC 
On the contrary, very little is known on the nuclear gluon density per nucleon, 
G {x,Q'^) (see [2] for a recent review). So far, the only constraint on G has 
been obtained from the scaling violation of the Fj" / ratio [3] measured 
by NMC [4]. The range of x explored in this experiment is 0.02-0.2, with a 
few GeV^-wide Q^-band around = x x 100 GeV^. However, the too large 
experimental uncertainties in these data do not allow for a precise estimate of 
the nuclear gluon distribution ratio, 

<(x,g2) = G^x,gyG^(x,g2). (i) 

The presently lack of knowledge on G therefore prevents reliable pQCD pre- 
dictions in high-energy nuclear reactions, such as p-Pb (a/s^ = 8.8 TeV) and 
Pb-Pb {^/s^ = 5.5 TeV) colhsions at the LHC. 

In addition to being a useful practical tool to predict hard processes in hadronic 
collisions, a precise knowledge of the gluon distribution is essential for a bet- 
ter understanding of evolution in QCD. In particular, an important theo- 
retical activity has recently focused on possible non-linear QCD evolution 
at small x and small Q^, where the gluon density in the nucleus becomes 
large, G (x, Q^) ~ \/as{Ql), and starts to saturate (Qs stands for "satura- 
tion scale") [5]. Measuring G in the vicinity of the saturation region, Q > Qs, 
would therefore provide useful tests of the saturation picture. 

A 

For these reasons, it is essential to determine G from the experimental data 
at RHIC and at the LHC with a similar accuracy to what has been achieved 
at HERA and Tevatron in the proton case. In p-A collisions, there are many 
ways to extract the nuclear gluon distribution G (x, Q"^) or the nuclear gluon 
distribution ratio R^. Let us discuss briefly various observables which look 
promising to achieve such a goal: 

• The prompt photon production channePI [6], explored in detail in this 
Letter, has a rich phenomenology, ranging from fixed-target experiments to 
the TevatroiH]. Next-to-leading order (NLO) pQCD calculation provides an 
impressive description of the world-data on > 5 GeV hadroproduction 
spectra [7] . As will be made more precise later, the mid-rapidity cross section 
is sensitive to parton distributions at x ~ x^ = 2^)^! \fs and ~ p^. At a 

^ In the following, "photons" always refer to prompt photons. In particular, we do 
not consider photons coming from hadron decays. 

^ At collider energies, the separation from 7r*^-decay photons becomes difficult with- 
out the use of isolation criteria. This issue will be discussed further in Section 3. 
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given Q^, other values of x can be probed by measuring photons at non-zero 

rapidity; 

• Jet production is a direct competitor to the photon channel which has 
the clear advantage of large counting rates at the LHC. The lowest p^^^ > 
30 GeV above which the jet signal becomes tractable is, nevertheless, rather 
large as compared to that of photons {p"' > 5 GcV), making both observ- 
ablcs complementary for the exploration of parton distributions. This is 
especially true for the typical Q^-range probed by those two mechanisms. 
However, since nuclear modifications of gluon densities prove more pro- 
nounced at low Q^, the photon process appears to be more appropriate to 
determine gluon distribution ratios; 

• Large-p^ dilepton production at small invariant mass, M < p^, has rather 
recently been put forward as a surrogate to photon production [8,9]. The 
main advantage of this channel is the absence of a large background. This 
channel is however suppressed by a and as compared to photons and 
jets, respectively, leading to much smaller counting rates; 

• Heavy boson {Z^, W^) production at the LHC has also been proposed 
to probe the shadowing of sea quarks (which are driven by gluons) at large 
scales ~ 10^ GeV^ [10, 11], even though nuclear effects are not expected 
to be too pronounced in this region. 

The kinematical (x, Q^)-domain probed with photons, jets, and heavy bosons 
(together with the range covered by NMC [4]) is displayed in Fig. 1, which 
highlights the complement ariness of these various probes. The dash-dotted 
line indicates the value of the saturation scale, Ql — 1 GeV^ (10^^/a;)°-^, 
obtained from the "geometrical scahng" fit to HERA data [12], scaled by 
^1/3 — 208^/^ in a Pb nucleus. As already pointed out, the photon channel 
covers in particular an important region at scales close to Qs where shadowing 
is strongest. 

The above list of observables is not meant to be comprehensive. In particular, 
more information could in principle be obtained from measurements of open 

heavy-flavoured mesons or large-p^ light-hadrons. As we shall see later, how- 
ever, the fragmentation process from partons to hadrons does not allow for 
the partonic kinematics to be determined. As a consequence, these processes 

A A 

rather probe moments of G rather than G itself. 

So far, none of the above channels has produced any constraint on G because 
of large error bars. At LHC, all these probes should be used to gain the 
most precise knowledge on parton distributions, and among them on G . First 
investigations fitting in such a program have been described in Ref. [13]. Below, 
a deeper analysis of the photon production channel is outlined. 
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Fig. 1. Typical (x = x^e~y,Q'^ = m^)-domain probed in p-Pb collisions at RHIC 
and LHC using photon, jet, and heavy boson production. At LHC, photons produced 
at = 5-100 GeV and \y\ < 3 are indicated, jets from = 30 GeV {\y\ < 4.5) to 
p^ = 300 GeV (|y| < 2), and heavy bosons at p^ < 100 GeV and \y\ < 3. At RHIC, 
photons produced at Pj^ = 3-6 GeV and y = 2-3 are selected. The kinematical range 
covered by the NMC experiment is also shown for comparison. The dash-dotted line 
indicates the saturation scale in a Pb-nucleus (see text for details). 

2 Inclusive photon production in p—A collisions 

2. 1 Perturbative dynamical processes 

At leading-order (LO), i.e. 0{aas), large-p^ prompt photon production pro- 
ceeds via 2 — i> 2 processes at the parton level. There are two channels in which 
the photon is produced directly in the partonic scattering: 

q{q) + G ^ 7 + q{q) (Compton), q + q •y + G (annihilation), 

yet the annihilation channel is less than a tenth of the Compton scattering 
below Xj_ < 0.1, hence the interest of photon production for probing gluon 
densities. At this order, the cross section for the direct channel is an integral 
over both the projectile and target parton momentum fractions, Xi and X2, of 
the elementary cross section multiplied by the parton distribution functions 
(PDF), generically fi^"\xi) for the projectile and fj^^^{x2) for the target, 
summed of the parton flavours i and j. Because of the 2 — 2 scattering 
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kinematics, the cross section reduces to a single integral along the hyperbola 

in the (xi, a;2)-plane. In order to understand the phenomenology discussed 
in the following, it is useful to have in mind some well-known basic trends 
deduced from the above equation. First, increasing Xj_ means exploring larger 

Xi,X2 (for instance the symmetric point. Xi = X2 ai y = reads Xi = X2 = x ^) . 
Next, going to positive y corresponds to probe larger Xi in the projectile 
(proton) and smaller X2 in the target (nucleus). 

In addition to the above direct process, a parton formed in a 2 — > 2 process 
can fragment into a coUinear photon, 

gi + g2 ^ gs + 94 ; ga ^ 7 + X , 

where q stands here for a generic parton. As thoroughly discussed in Ref. [14], 
although the 2 — > 2 scattering process is formally O (ct^), the parton-to-photon 
fragmentation mechanism O (a/ag) makes this channel to contribute also at 
leading-order. The cross section in this channel has the same ingredients as 
the direct process, supplemented by the fragmentation process. The latter is 
encoded in a fragmentation function Dg.^^^{z), where the momentum-fraction 
z of the parent-parton carried away by the photon has to be integrated from 
x^coshi/ to 1. Consequently, the production of a fragmentation photon at 
a given and y probes parton distributions at larger momentum-fractions 
{xi/z,X2/z) than at {xi,X2) in the direct channel. This can be read off the 
hyperbola equation, after the change x^/2 x^/{2z) is implemented to take 
into account the effect of fragmentation in the kinematics. 

Exploration of nuclear modifications to parton densities in such a framework 
simply consists in the replacement, say for the target, of the proton densities 
by the nuclear parton distribution functions (nPDF): 

fr{x2) ^ ftix2). 

Since the cross section is actually known at NLO, 0{aal), the analysis will 
be carried out at this level of accuracy. 

2.2 Nuclear production ratio at the LHC 

The inclusive photon - that is summing direct and fragmentation channels - 
p^ spectra has been computed in p-p and p-A collisions at 8.8 TcV (LHC 
conditions) using the work of Ref. [14] . In p-p scattering, the NLO CTEQ6M 
parton densities have been used in the calculation [1] , and the fragmentation 
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functions of quarks and gluons into photons were taken from the NLO fit of 
e"'"e~ data carried out in Ref. [15]. Regarding nuclear PDFs, is obtained 
from the average over the proton and neutron distributions in a nucleus with 
atomic mass A and atomic number Z, ff" = Z ff^^ + {A — Z) fi^^. Several 
parametrizations of the ratio of the proton distribution in a nucleus, ff^^, over 
that of a "free" proton, 

<(^,Q')=/f^^x,Q2)//f(x,Q^), (2) 

are available for each parton flavour i: EKS [16], HKM [17], nDS [18]; yet only 
the latter group has performed a global data analysis at NLO accuracy, used 
in the present LetteiZl- An alternative fit, for which the depletion of gluon 
shadowing is somehow arbitrarily enhanced at small x as compared to nDS, 
is proposed in [18] and also considered for comparison (labelled nDSg in the 
following) . 

In order to quantify nuclear effects, the nuclear production ratio 



is determined at LHC in the = 5-100 GeV range and at rapiditjQ 1/ = 0. 
It is plotted in Fig. 2 as a function of . i^^p^^ smoothly increases from lower 
to higher x^. At x^^ ~ 10~^, the suppression is only roughly 10% (-R^pj^ — 
0.92) when using nDS parton densities, but proves larger {R^p^^ — 0.85) if the 
stronger nDSg gluon shadowing is assumed in the calculation. At large Xj_ 
(recall this would correspond to large hence to large Q^), there is basically 
no attenuation of the photon yield in p-Pb with respect to p-p collisions. The 
trend of the inclusive photon suppression versus x^ can be simply understood: 
as x^ gets larger, the partons probed in the nucleus carry a higher momentum- 
fraction X2 for which the gluon nuclear density ratio R is less suppressed (i.e. 
closer to 1). Another reason comes from the Q^-dependence of the shadowing 
process. Indeed, large-Xj^ photons are produced perturbatively from highly- 
virtual partons = O (jP^^ (at fixed energy, x^ and are correlated) that 
are less affected by nuclear shadowing corrections, which die out at asymptotic 

Similarly, this calculation is performed for inclusive photons produced at ra- 
pidity y = 2.5. However, since the results are rather close to what is found 



^ For completeness, we checked that the mam conclusion of this study does not 
depend much on which parametrization is used in the calculation. 
^ In a p-Ph at ^/s^ = 8.8 TeV at LHC, this corresponds to y^^^ ^ 0.5 in the 
laboratory frame. However, we checked that the results are practically unaffected 
by the small rapidity-shift from the centre-of-mass to the laboratory frame. 
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in the isolated photon channel discussed later, photon production at forward 
rapidities is discussed in Section 3. 
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Fig. 2. Nuclear production ratio ^pp,^ of inclusive photon production at y = in 



p-Ph collisions at 



8.8 TeV. 



Let us now briefly discuss to which extent such a suppression could be seen 
at the LHC At mid-rapidity, the differential absolute cross section in p-Ph 
collisions is roughly 



daP- 



-Pb 



dy dp^ 
dcjP-P^ 



dy dp^ 



=10 GeV 



p . =100 GeV 



1.4 10^ pb/GeV, 



0.9 10^ pb/GeV. 



(4) 



This would correspond to a yearly rate of A/" ~ 10^ and 10^ per GeV-bin and 
per unit rapidity, at p^^ = 10 and 100 GeV respectively, with the expected 
integrated luminosity in p-Pb collisions at the LHC £ = 1.4 10^^ cm~^ (tak- 
ing 1 year = 10^ s) [13]. Neglecting the statistical errors in the p-p run taken 
at high luminosity, the statistical uncertainty on the nuclear production ratio 
is therefore negligible at p^ = 10 GeV and roughly SR^^^ ~ A/"^^/^ ~ 3% 
in a 1 GeV-bin at p^ = 100 GeV. Perhaps with the exception of the mod- 
erate suppression when using nDS for mid-rapidity photons, the fact that 



R in all other cases indicates that the expected attenuation 



pPb 

of inclusive photon at small could be observed experimentally, as far as 
statistical signiflcance is concerned. Of course, this rough estimate does not 
include any consideration on detector acceptanc^ and efficiency, nor possi- 
ble systematic errors, which should be studied in detail in the experimental 
analysis. Nevertheless, we find it encouraging that, despite the much smaller 
photon rate than that of jets, the cross sections at the LHC are large enough 



^ Note for instance the rather limited coverage of the PHOS electromagnetic 
calorimeter in the ALICE experiment, Ay x A(j) = 0.26 x 1.7. 
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to measure nuclear effects in this channel up to ~ 10^ GeV {Xj_ ~ 2 10 ). 
Let us also mention that a similar accuracy is expected for photons produced 
at forward rapidity {y = 2.5), which could be measured e.g. by the CMS 
experiment [19]. 



3 Extracting gluon nuclear distributions 

3. 1 Isolated photons and the relationship with gluon distributions 

We now turn to the study of isolated photon cross section. In this case, there 
is a very simple relationship between R^^{x^) and the parton distribution 
ratios, Eq. (2), R^{x,Q'^). This relationship stems from the slow variation of 
the latter ratios as compared to those of the parton distributions themselves, 
ff'{x,Q^). It is approximate, but the deviation from the exact result may 
well be hidden by other sources of uncertainty, in which circumstances such a 
direct extration is most valuable. As a matter of fact, the approximation can 
be gauged a posteriori: the weaker the x^-dependence of the measured ratio, 
-^pa(^x)) the better the approximation. 

The LO cross section is a sum of convolutions of projectile and target par- 
ton densities with partonic cross sections. Writing xi = x_[_e'^/{2v) the LO 
Compton channel cross section reads: 

dV ^ a as ta..g / x^e'^ \ 1 / 1 \ 

dyd?p^ {x^s/2)^ \ 2v J [2{l-v))3\ 1-v) 

where F{x) = F^{x)/x. For proton-nucleus interaction, the per-nucleon cross 
section has the same expression with the changes 

f^'\x) ^ F^{x) = b!1^{x)F\x), G'"'(x) ^ G^{x) = rI{x)G\x), 

where we recall that F^ and are understood as per-nucleon distribution^^. 

The integrand in (5) is strongly suppressed at the end-points. This is a result 
of the competition between F'^^°^ and G""^^ (as well as and F'"''^) that 
are suppressed at large values of their respective arguments. Since the PDF 

^ There is also a change Vmax — > 1 — Xj_e~y/{2A), which is of no practical impor- 
tance. 
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ratios, R^^ and R^, have a slow variation as compared to the F x G product, 
it can be assumed to be constant and put out of the w-integral. The typical 
X at which R is probed is easy to estimate at small- and for \y\ not too 
large, where the integrand is driven by the small-x dependence of F's and G's. 
Considering F{x) ~ Ax~"' and G{x) ~ Bx~^ at small x, quantities such as 
— v)^ show up for F X G products. If the difference between a and b is 
disregarded, F x G is maximal at f = 1/2 (a, 6 > 0). This is not the whole 
story, since the f-dependence beside that oi F x G shifts the maximum of the 
integrand away from v = 1/2. However, this is partly counterbalanced by the 
fact that a < b (around x — 10~^, a 1.3 and b ^ 1.7). In the following the 
simple replacement 




will be studied. 

Making this replacement in Eq. (5) leads to a very simple relationship between 

A A 

R^j^ and nuclear distribution ratios: R^^ is an average of R^ and R^ weighted 
by y-dependent coefficients that come from the relative importance of the two 
terms in Eq. (5) evaluated for p-p. Rather than quoting the full result, three 
simple cases of interest will be discussed. At y = 0, thanks to the symmetry 
between the two terms in Eq. (5) for p-p, the weights are 0.5, hence 

R^,ix,,y = 0) ^ 0.5 <Ja;J + 0.5 R^^ixJ. (6) 



At large y > yo (in practice yQ ~ 2.5), i.e. at somewhat larger xi (notice 
that the above reasoning at small x is not adequate to properly describe this), 
the extinction of gluons leads to a suppression of the second term in Eq. (5), 
leading to 

R^A^,,y>yo)^Ki^^e-y). (7) 

Conversely, at y < —yo, 

RAi^.,y<-yo)^K^{x^ey). (8) 



3.2 Results at LHC and RHIC 



In order to check the accuracy of the assumptions of the former Section, the 
production of isolated photons is computed in p-Pb collisions at ^/s^ — 
8.8 TeV in pQCD at NLO, using nDSg nuclear parton densities. The isolation 
criterion was chosen as follows: a photon is said to be isolated if the total 
transverse partonic energy inside a cone of radius R = \J Ay^ + = 0.4 
around the photon is less than 10% of its transverse momentum: Ej_^^^^{R ~ 
0.4) < 0.1 p.. 
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The nuclear production ratio at ?/ = and y = 2.5 is plotted as a function of 
Xj^e~^ in Fig. 3 and Fig. 4. The -dependence of R^^^ at |/ = is similar to 
what is observed in the inclusive channel, yet the suppression of the isolated 
photon yield, especially at low , turns out to be slightly more pronounced. 
At mid-rapidity, the production ratio above a;^ ^ 3 10^^ is consistent with 
1, indicating small nuclear effects in the parton distributions. The forward- 
rapidity photon suppression is larger, at a given Xj_, than at mid-rapiditjQ 
and does not vanish even in the largest x^-bin, x^e~^ ~ 4 10~^. 
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Fig. 3. Nuclear production ratio -Rpp^ 
of isolated photon production at y = 



in p-Ph collisions at 



TeV. 



Fig. 4. Same as Fig. 3 for photons with 
rapidity y = 2.5. 



Let us now compare the isolated photon suppression, R^^, with the analytic 
approximations R^'''"'' , Eqs. (6) to (8), in terms of parton densities ratios. 

In Fig. 3, R^^q''{x^) = 0.5 R^{x^) + -R^^(x^) is plotted as a dotted line. 
The agreement between this approximation and the exact suppression ratio 
(solid line) is very good, of the order of one-percent accuracy on the full 



x^ range, as can be seen from the ratio r^=^ = {R^j^ — R^y^^"^) / R^^^ plotted 



as a dash-dotted line. Therefore, it turns out that neither the annihilation 
process nor next-to-leading order corrections, both neglected in Section 3.1, 
spoil dramatically the quahty of the analytic approximation. Consequently, the 
nuclear production ratio of isolated photons can serve as a reliable probe to 
measure small-x shadowing of parton densities in large nuclei. The production 
oi y = 2.5 isolated photons shows a similar trend. At large rapidity, R^^ 



can be approximated with R, 



-.approx, - _ A. 

j>yo y^±J ~ ^G^^J 

line), together with r^'^^ = {Rp^- R7>17)/ (dash-dotted line), showing an 



e-y) plotted in Fig. 4 (dotted 



If, on the contrary, one compares the suppression at a given x' = x^e ^, the 
forward-rapidity photons are less suppressed than at y = 0. Indeed, the larger the 



rapidity, the more virtual — "'^^ ""''^ 
and therefore the lesser the nuclear shadowing. 



x'^ sjAi the partons are probed in the nucleus, 
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agreement better than 5%. At both rapidities, a similar behaviour is observed: 
R^j^ is shghtly above gx small x and then becomes closer to /^"'''"'™ 

large Xj_. The fact that R^j^ > partly comes from the extra gluon 

radiated at NLO which requires the momentum-fraction carried by the target 
parton to be larger than it would at LO accuracy. 

We now present the production ratio of isolated photons in d-Au (over p-p) 
collisions at RHIC top-energy, y^s^ = 200 GeV. The calculation has been 
carried out in the forward-rapidity region [y = 3) in the p^^ = 3-6 GeV range, 
corresponding to Xj_e'~y — 1.5 10~^-3 10"^. Note that measurements at such 
large rapidities, y ^ 3, could be envisaged both with the PHENIX and the 
STAR detectors [20]. 

Before we go on, let us stress that a reliable pQCD calculation can be seen 
as doubtful for two reasons. First of all, in that small p^^-domain (basically 
dictated by kinematic constraints), ag is not too small and PDFs and FFs are 
not well constrained at such low virtualities. Moreover, the present fixed-order 
calculation might be somehow affected by large threshold and recoil resum- 
mation corrections at the edge of phase-space, that is at large p^ [21]. That 
said, let us explore how good the approximated expressions for i?^^^ might be 
at RHIC within the current approach. The ratio R^j^^, plotted in Fig. 5 (solid 
hue) together with the estimate R'^yJ^ixj^) — R^,{x^e^y) (dashed), does not 
depend on x^ in the limited coverage assumed in the calculation. The gluon 
distribution ratio, i?^ ~ 0.8, is similar to what is found at LHC and lie some- 
what above i?^^^ . The disagreement between i?^^^ and R^J>y^ is actually due to 
a trivial isospin effect in the deuteron projectile: the presence of the neutron, 
less efficient than a proton to produce prompt photons, leads to an additional 
suppression. This effect can however easily be corrected for, e.g. by multiplying 
^dAu t»y the inclusive production ratio in d-p over p-p collision. This isospin- 
corrected ratio (labelled R^ * dp/pp and displayed as a dotted line) proves 
close to the measured R^^^ ratio, rjg^^ = (i?,^^ - i?^> J * dp/pp) /R^^^ ~ 8%, 
yet the agreement is somewhat less than at LHC energy. 



3.3 Extraction without a reference to p-p 

The direct connection between production ratios and nuclear distribution ra- 
tios is of course interesting. It shows the need for having the cross section in 

p-p and p-A in the same kinematical conditions. At LHC, this necessitates an 
extrapolation in energy (the nominal center of mass energy is 14 TeV in p-p, 
and 8.8 TeV in p-Pb) and a comparison of event rates at different rapidities in 
the detector (in p-Pb the center of mass frame moves at a rapidity 0.47 in the 
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Fig. 6. Ratio of forward-over-back- 
ward cr{+y)/a{—y) inclusive photon 
production at y = 2.5 in p-Ph colli- 



sions at 



8.8 TeV. 



detector frame). This gives extra sources of uncertainties in the extractio 



There is a way out of this difficulty which can be obtained by comparing in 
the same experiment the p-A cross section at forward {+y) and backward 
{—y) rapidity. Thanks to the symmetry of the p-p cross section in the change 
y —y, the ratio 



y) _RJx,e-y) 



the latter relationship coming from the approximations at \y\ 
cussed previously. At large y, there is thus a way of obtaining at x^e~y 



> 2.5 dis- 

A 

from the knowledge of i? at x , e"^^. As an illustration, the forward-over- 



backward ratio in p^Pb coflisions is plotted in Fig. 6 together with the esti 

t 

as long as is not too large (x^e~^ < 10"'') 



y) / R^^{x ^6^). This approximation turns out to be good (~ 5%) 

y < 10~3^ 



At LHC, thepj^ range 5-100 GeV corresponds to Xj_e~^y ^ 0.01-0.03 aty = 2.5, 
that is precisely the NMC x-range (see Fig. 1). Making the Q^-evolution of 
the NMC ratio to = p2 (^^^^ 

in a larger nucleus) leads to the following 



^ In that sense, performing measurements in d-d and d-Pb collisions (and to a lesser 
extent in p-d and p-Ph) would have the double advantage to get rid of possible 
isospin effects when comparing nuclei with very different Z/A ratios and to avoid 
the boost of the center-of-mass frame with respect to that of the laboratory. 
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estimate for 



K(-.--',pI) - ^^^f X flr(-.»".!'!) i« > 2.5) 



4 Summary 



There is a major importance to determine precisely the nuclear gluon distri- 
butions at small x, whether it is to perform accurate pQCD predictions or 
to study QCD evolution equations. We explored in this Letter how prompt 
photon production in p-Ph collisions at LHC and d-Au collisions at RHIC 
can help to achieve this goal. 

The nuclear production ratio of isolated photons in p-A collisions can be sim- 
ply approximated as a linear combination of gluon distributions and structure 

A A 

functions in the nucleus A over those in a proton, and R^ . Performing 
the calculation in pQCD at NLO we checked that such an approximation is 
correct up to a few-percent accuracy, making this observable an ideal tool to 
measure gluon shadowing at RHIC and at LHC. We also point out that the 
ratio R^{x^ e~'^)/R^ {x^ e^) can be determined through the forward-over- 
backward photon production ratio in p-A collisions, without the need of any 
p-p reference data at the same energy. 

Although a particular set of the nuclear parton densities is used here as an 
example, the method and results obtained in this study are quite general. 
Present fits and models of nPDFs are not constrained in the small (x, Q"^)- 
domain, with a factor of 2 spread at x ~ 10^^ [2]; the nDSg set being one 
of the fits with the weakest shadowing. Future prompt photon data, therefore 
with a suppression potentially larger than what is determined in this Letter, 
should be able to properly discriminate among the parameterizations currently 
available. 
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